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Abstract: On directed Barabasi-Albert networks with two and seven neigh- 
bours selected by each added site, the Ising model was seen not to show a spon- 
taneous magnetisation. Instead, the decay time for flipping of the magnetisation 
followed an Arrhenius law for Metropolis and Glauber algorithms, but for Wolff 
cluster flipping the magnetisation decayed exponentially with time. On these net- 
works the Majority-vote model with noise is now studied through Monte Carlo 
simulations. However, in this model, the order-disorder phase transition of the or- 
der parameter is well defined in this system. We calculate the value of the critical 
noise parameter q c for several values of connectivity z of the directed Barabasi- 
Albert network. The critical exponentes (3/v, 7/V and 1/v were calculated for 
several values of z. 
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Introduction 

It has been argued that nonequilibrium stochastic spin systems on reg- 
ular square lattice with up-down symmetry fall in the universality class of 
the equilibrium Ising model [1]. This conjecture was found in several models 
that do not obey detailed balance [21 El El- Campos et al. |H] investigated the 
majority- vote model on small-world network by rewiring the two dimensional 
square lattice. These small-world networks, aside from presenting quenched 
disorder, also posses long-range interactions. They found that the critical 
exponents j/u and (3/v are different from the Ising model and depend on the 
rewiring probability. However, it was not evident that the exponent change 
was due to the disordered nature of the network or due to the presence of 
long-range interactions. Lima et al. jH] studied the majority-vote model on 
Voronoi-Delaunay random lattices with periodic boundary conditions. These 
lattices posses natural quenched disorder in their conecctions. They showed 
that presence of quenched connectivity disorder is enough to alter the ex- 
ponents (3/v and 7/1/ from the pure model and therefore that is a relevant 
term to such non-equilibrium phase-transition. Sumour and Shabat [3 |S] 
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investigated Ising models on directed Barabasi-Albert networks jHj with the 
usual Glauber dynamics. No spontaneous magnetisation was found, in con- 
trast to the case of undirected Barabasi-Albert networks [TO! ITT| IT2~] where 
a spontaneous magnetisation was found lower a critical temperature which 
increases logarithmically with system size. More recently, Lima and Stauffer 
ilSj simulated directed square, cubic and hypercubic lattices in two to five 
dimensions with heat bath dynamics in order to separate the network ef- 
fects from the effects of directedness. They also compared different spin flip 
algorithms, including cluster flips jT3], for Ising-Barabasi- Albert networks. 
They found a freezing- in of the magnetisation similar to [3 |K| , following an 
Arrhenius law at least in low dimensions. This lack of a spontaneous mag- 
netisation (in the usual sense) is consistent with the fact that if on a directed 
lattice a spin Sj influences spin Si, then spin Si in turn does not influence 
Sj, and there may be no well-defined total energy. Thus, they show that for 
the same scale-free networks, different algorithms give different results. Now 
we study the Majority-vote model on directed Barabasi-Albert network and 
different from the Ising model, the order- disorder phase transition of order 
parameter well it is defined in this system. We calculate the (3/v, j/u, and 
1/u exponents and these are different from the Ising model and depend on 
the values of connectivity z of the directed Barabasi-Albert network. 



Model and Simulaton 

We consider the majority-vote model, on directed Barabasi-Albert Net- 
works, defined [T5J EH EH UZj by a set of "voters" or spins variables o taking 
the values +1 or —1, situated on every site of a directed Barabasi-Albert 
Networks with N sites, and evolving in time by single spin-flip like dynamics 
with a probability Wi given by 



Wi( a ) = g 



I -(l-2q)a i S[J2a r+5 

■8=1 



(1) 



where S(x) is the sign ±1 of x if x ^ 0, S(x) — if x — 0, and the sum 
runs over all nearest neighbors of <7j. In this network, each new site added to 
the network selects z already existing sites as neighbours influencing it; the 
newly added spin does not influence these neighbours. The control parameter 
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q plays the role of the temperature in equilibrium systems and measures the 
probability of aligning antiparalle to the majority of neighbors. 

To study the critical behavior of the model we define the variable m = 
YliLi&i/N. In particular , we were interested in the magnetisation, suscep- 
tibility and the reduced fourth-order cumulant: 

M(q) = [< \m\ >] (2) 

X (q) = N[<m 2 >-< \m\ > 2 ] av , (3) 



U{q) = 



< m 4 > 



(4) 



3 < m 2 > 2 

where < ... > stands for a thermodynamics average and [...] a „ square brackets 
for a averages over the 20 realizations. 

These quantities are functions of the noise parameter q and obey the 
finite-size scaling relations 

M = N-M"f m (x)[l + ...], (5) 
X = N^f x (x)[l + ...], (6) 

^ = N 1 '»f u (x)[l + ...], (7) 

where u, (3, and 7 are the usual critical exponents, fi(x) are the finite size 
scaling functions with 

x=(q- q c )N^ (8) 

being the scaling variable, and the brackets [1 + ...] indicate corretions-to- 
scaling terms. Therefore, from the size dependence of M and \ we obtained 
the exponents (3/v and 7/^, respectively. The maximum value of susceptibil- 
ity also scales as N 1 ^ . Moreover, the value of q for which \ has a maximum, 
qXmax — q c (N), is expected to scale with the system size as 

q c (N) =q c + bN- 1 /", (9) 

were the constant b is close to unity. Therefore, the relations (7) and (9) are 
used to determine the exponente I /v. We have checked also if the calculated 
exponents satisfy the hyperscaling hypothesis 

2/3/u + 7 /u = D eff (10) 
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in order to get the effective dimensionality, D e ff, f° r various values of z. 

We have performed Monte Carlo simulation on directed Barabasi-Albert 
networks with various values of connectivity z. For a given z, we used systems 
of size N = 1000, 2000, 4000, 8000, and 16000. We waited 10000 Monte Carlo 
steps (MCS) to make the system reach the steady state, and the time averages 
were estimated from the next 10000 MCS. In our simulations, one MCS is 
accomplished after all the N spins are updated. For all sets of parameters, 
we have generated 20 distinct networks, and have simulated 20 independent 
runs for each distinct network. 

Results and Discussion 

In Fig. 1 we show the dependence of the magnetisation M and the sus- 
ceptiblity x on the noise parameter, obtained from simulations on directed 
Barabasi-Albert network with 16000 sites and several values of conectivity z. 
In the part (a) each curve for M, for a given value of N and z, suggests that 
there is a phase transition from an ordered state to a disordered state. The 
phase transition occurs at a value of the critical noise parameter q c , which 
is an increasing function the conectivity z of the directed Barabasi-Albert 
network. In the part (b) we show the corresponding behavior of the suscep- 
tibility Xi the value of q where \ has a maximum is here identified as q c . In 
Fig. 2 we plot the Binder's fourth-order cumulant for different values of N 
and two different values of z. The critical noise parameter q c , for a given 
value of z, is estimated as the point where the curves for different system 
sizes N intercept each other. In Fig 3. the phase diagram is shown as the 
dependence of the critical noise parameter q c on connectivity z obtained from 
the data of Fig. 2. 

The phase diagram of the majority-vote model on directed Barabasi- 
Albert network shows that for a given network (fixed z ) the system becomes 
ordered for q < q c , whereas it has zero magnetisation for q > q c . We notice 
that the increase of q c as a function of z is slower than the one than in [T7| . 
In the Fig. 4 we plot the dependence of the magnetisation at q = q c with the 
system size. The slopes of curves correspond to the exponent ratio (3/v of 
according to Eq. (5). The results show that the exponent ratio (3/v decreases 
when z increases, see Table I. 

In Fig. 5 we display the scalings for susceptibility at x(lc(N)) (circle) 
for its maximum amplitude Xtv^S an d x{Qc{N)) (square) obtained from the 
Binder's cumulant versus N for connectivity z = 8. The exponents ratios 
j/v are obtained from the slopes of the straight lines. For almost all the 
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values of z, the exponents j/u of the two estimates disagree (Table I). An 
increased z means a tendency to increase the exponent ratio 7/V, see Table 
I, so that they disagree with the results of Luiz et al [T7] , where the values 
of exponents ratio 7/V are almost all equal and with a slight tendency to 
decrease. Therefore we cannot use the Eq.(9), for fixed z, obtain the critical 
exponent 1/V. In the Fig. 6 we show the critical behavior of /3/V and 7/V as 
a function of connectivity z. 

To obtain the critical exponent we calculated numerically U (q) = 
dU(q)/dq at the critical point for each values of iV at connectivity fixed z. 
The results are bad agreement with the scaling relation (7). Then, also we 
cannot calculate the exponents 1/V, through this relation. Therefore we do 
not obtain to get the values of the exponents 1/V for each connectivity z 

The Table I resumes the values of q c , the exponents /3/v, 7/1/, and the 
effective dimensionality of systems. For all values of z the value -D e // = 1, 
which has been obtained from the Eq. (9), therefore when z increases, (3/v 
decreases and 7/1/ increases, thus providing the value of -D e // = 1 (along 
with errors). Therefore, the directed Barabasi- Albert network has the same 
effective dimensionality as Erdos-Renyi's random graphs |T7j . J. M. Oliveira 
[To] showed which majority-vote model defined on regular lattice has critical 
exponents that fall into the same class of universality as the corresponding 
equilibrium Ising model. Campos et al [S] investigated the critical behavior of 
the majority- vote on small- world networks by rewiring the two-dimensional 
square lattice, Luiz et al ^7] studied this model on Erdos-Renyi's random 
graphs, and Lima et al [0] also studied this model on Voronoi-Delaunay lat- 
tice. The results obtained by these authors show that the critical exponents 
of majority-vote model belong to different universality classes. 

Finally, we remark that our MC results obtained on directed Barabasi- 
Albert network and undirected (in preparation) majority-vote model show 
that critical exponents are different from the results of for regular lattice 
and of Luiz et al (T7j for Erdos-Renyi's random graphs . 

Conclusion 

In conclusion, we have presented a very simple nonequilibrium model 
on directed Barabasi-Albert network [7||Hj. Different from the Ising model, 
in these networks, the Majority- vote model presents a second-order phase 
transition which occurs in model with connectivity z > 1. The exponents 
obtained are differentfrom the other models. Nevertheless, our Monte Carlo 
simulations have demonstrated that the effective dimensionality -D e // equals 
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D eff 
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n a^kacw 




U.UU'i^O J 




l.UlOl £? J 


3 


0.431(4) 


0.444(1) 


0.111(2) 


0.904(12) 


0.999(2) 


4 


0.431(3) 


0.447(1) 


0.104(2) 


0.888(9) 


0.998(3) 


6 


0.438(2) 


0.435(2) 


0.134(5) 


0.861(3) 


1.008(6) 


8 


0.444(5) 


0.431(1) 


0.138(2) 


0.851(5) 


1.000(2) 


10 


0.446(3) 


0.421(2) 


0.158(3) 


0.834(7) 


1.000(5) 


20 


0.458(4) 


0.412(1) 


0.178(2) 


0.795(11) 


1.002(2) 


50 


0.467(2) 


0.375(4) 


0.249(7) 


0.735(17) 


0.999(11) 


100 


0.474(3) 


0.363(4) 


0.269(5) 


0.674(23) 


0.999(9) 



Table 1: The critical noise q c , the critical exponents, and the effective di- 
mensionality D e ff , for directed Barabasi- Albert network with connectivity 
z . 

units, for all values of z, that are agree with the results de Luiz et al [T7] . 
However, when z grows, the exponents at the critical point q c , f3/v obtained 
by Binder's cumulant decrease and the exponents j/u grow, satisfying the 
hyperscaling relation with D e ff = 1. 
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Figure 1: Magnetisation and susceptibility as a function of the noise param- 
eter q, for N = 16000 sites. From left to ri,ht, z = 2, 3, 4, 6, 8, 10, 20, 50, 
and 100 . 
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Figure 2: Binder's fourt-order cumulant as a function of q. In part (a) we 
have z = 3 and part (b) z = 50. 
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Figure 3: The phase diagram, showing the dependence of critical noise pa- 
rameter q c on connectivity z. 
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Figure 6: Critical behavior the (3/v and 7/V exponents as a function of 
connectivity z. 
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